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We present results on time-dependent CP asymmetries in neutral B decays to several CP eigen- 
states. The measurements use a data sample of about 227 x 10 6 T(4S) — > BB decays collected by 
the BABAR detector at the PEP-II asymmetric-energy B Factory at SLAC. The amplitude of the 
CP asymmetry, sin2/3 in the Standard Model, is derived from decay-time distributions from events 
in which one neutral B meson is fully reconstructed in a final state containing a charmonium meson 
and the other B meson is determined to be either a B° or B° from its decay products. We measure 
sin2/3 = 0.722 ± 0.040(stat) ± 0.023(syst) in agreement with the Standard Model expectation. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



Charge-parity (CP) violation in the B meson system 
has been established by the BABAR 1] and Belle 
collaborations. The Standard Model of electroweak in- 
teractions describes CP violation as a consequence of 
an irreducible phase in the three-generation Cabibbo- 
Kobayashi-Maskawa (CKM) quark-mixing matrix Q . In 
this framework, measurements of CP asymmetries in 
the proper-time distribution of neutral B decays to CP 
eigenstates containing a charmonium and K° meson pro- 
vide a direct measurement of sin2/3 Q. The angle (5 is 
arg[-V cd V* h /V td V t l], where Vy are CKM matrix ele- 
ments. 

In this Letter we report on an updated measurement 
of sin2/3 in (227 ± 2) x 10 6 BB decays using B° decays to 
the final states J/ipK° s , J/ipK° ip(2S)K° s , XdK° s , rj c K° s , 
and J/ipK* a {K*° -> K°tt°) |J. The BABAR detector 
and the measurement technique are described in detail 
in Refs. Q and Q, respectively. Changes in the analysis 
with respect to the previously published result include 
140 x 10 6 more BB events, an improved event recon- 
struction applied to all of the data, a new flavor-tagging 
algorithm, and fewer assumptions about the CP proper- 
ties of background events. 

The proper-time distribution of B meson decays to a 
CP eigenstatef can be expressed in terms of a complex 
parameter A [8j , which depends on both the B°-B° oscil- 
lation amplitude and the decay amplitudes for B° — > / 
and B° — >■ /. The decay rate f+(f-) when the other B 
meson B tag decays as a B° (B°) is given by 



f±(Ai) = 



,-|At|/T B o 

4t r o 



2 Zm A . . . . . 
1 ± - rr^r sin (Am d At) 



1 + 



cos (Am^At) 



, (1) 



for a B from a T(AS) — ► B°B° decay, where At is the 
difference between the proper decay times of the recon- 
structed B meson B rec and -Btag, t b q is the B° lifetime, 
and Arrid is the B°-B° oscillation frequency. The decay 
width difference Ar between the B° mass eigenstates is 
assumed to be zero. The sine term is due to the interfer- 
ence between direct decay and decay after a net B°-B° 
oscillation. A non-zero cosine term arises from the in- 
terference between decay amplitudes with different weak 
and strong phases (direct CP violation) or from CP vio- 
lation in B°-B° mixing. 



In the Standard Model, CP violation in mixing is neg- 
ligible, as is direct CP violation for b — > ccs decays that 
contain a charmonium meson With these assump- 
tions A = r]fe~ 2113 , where rjf is the CP eigenvalue of final 
state /. Thus, the time-dependent CP asymmetry is 



Acp(At) 



f+-f- 
f++f_ 



-77/ sin2/3 sin (Am d At), (2) 



with r] f = -1 for J/il>K°, ip(2S)K°, XaK° s , and rj c K° s , 
and +1 for J/tpK®. Due to the presence of even (L=0, 
2) and odd (L=l) orbital angular momenta in the B — » 
J/i/jK*° final state, there can be CP-even and CP-odd 
contributions to the decay rate. When the angular infor- 
mation in the decay is ignored, the measured CP asym- 
metry in J/tp K*° is reduced by a factor 1 1 — 2R± | , where 
R± is the fraction of the L=l contribution. We have 
measured R± = 0.230 ±0.015 ±0.004 Q, which gives an 
effective 77/ = 0.51 ± 0.04, after acceptance corrections. 

In addition to the CP modes described above, we 
utilize a large sample (Bfl av ) of B° decays to the fla- 
vor eigenstates D^*^h + {h + = tt + , p + , and a~l) and 
J/ip K*°(K*° — > A'+7r~) for calibrating our flavor tag- 
ging and At resolution. Validation studies are per- 
formed with a control sample of B + mesons decaying 
to the final states J/ipK^+, il>(2S)K + , XdK + , and 
i] c K + . The event selection and candidate reconstruction 
are unchanged from those described in Refs. Q, 0, liot - 
except that only the 77 c — > K®K + ir~ channel is used 
in the B° -> r) c K° s and B ± -> r7 CJ Fr ± modes (2.91 < 
m K o K+7r - < 3.05GeV/c 2 ). 

The time interval At between the two B decays is cal- 
culated from the measured separation Az between the 
decay vertices of B rec and B tag along the collision (z) 
axis 0- We find the z position of the B rec vertex from 
its charged tracks. The B tag decay vertex is determined 
by fitting tracks not belonging to the B roc candidate to 
a common vertex, employing constraints from the beam 
spot location and the B lcc momentum 0- We accept 
events with a calculated At uncertainty of less than 2.5 ps 
and |Ai| < 20 ps. The fraction of events satisfying these 
requirements is 95%. The r.m.s. At resolution is 1.1 ps 
for the 99.7% of these events that exclude outliers. 

We use multivariate algorithms to identify signatures 
of B decays that determine ( "tag" ) the flavor at decay of 
the Btag to be either a B° or B°. Primary leptons from 
scmileptonic B decays are selected from identified elec- 



5 



TABLE I: Efficiencies t%, average mistag fractions Wi, mistag 
fraction differences Awi = Wi(B°) — Wi(B°), and Q extracted 
for each tagging category i from the Bfl av sample. 



Category 


e(%) 


w (%) 


Aw (%) 


Q (%) 


Lepton 


8.6 ± 0.1 


3.2 ± 0.4 


-0.2 ± 0.8 


7.5 ± 0.2 


Kaon I 


10.9 ± 0.1 


4.6 ± 0.5 


-0.7 ± 0.9 


9.0 ± 0.2 


Kaon II 


17.1 ± 0.1 


15.6 ± 0.5 


-0.7 ± 0.8 


8.1 ± 0.2 


Kaon-Pion 


13.7 ± 0.1 


23.7 ± 0.6 


-0.4 ± 1.0 


3.8 ± 0.2 


Pion 


14.5 ± 0.1 


33.0 ± 0.6 


5.1 ± 1.0 


1.7 ± 0.1 


Other 


10.0 ± 0.1 


41.1 ± 0.8 


2.4 ± 1.2 


0.3 ± 0.1 


All 


74.9 ± 0.2 






30.5 ± 0.4 



trons and muons as well as isolated energetic tracks. The 
charges of identified kaon candidates define a kaon tag. 
Soft pions from D* + decays are selected on the basis of 
their momentum and direction with respect to the thrust 
axis of £>tag- These algorithms are combined to account 
for correlations among different sources of flavor informa- 
tion and to provide an estimate of the mistag probability 
for each event. These algorithms have been improved 
relative to Ref. Q] with the addition of information from 
low-momentum electrons, A — > p7r decays, and additional 
correlations among identified kaon candidates. 

Each event is assigned to one of six tagging categories 
if the estimated mistag probability is less than 45%. The 
Lepton category contains events with an identified lep- 
ton; the remaining events are divided into the Kaon I, 
Kaon II, Kaon-Pion, Pion, or Other categories based on 
the estimated mistag probability. This new definition of 
tagging categories improves the overall performance of 
the tagging algorithm, while largely preserving the sepa- 
ration of events with differing sources of tagging informa- 
tion. For each category (i), the tagging efficiency £$ and 
fraction Wi of events having the wrong tag assignment 
are measured from data (Table |]J . The effective tagging 
efficiency Q = Y2i — 2u^) 2 improves by about 5% 
(relative) over the algorithm used in Ref. |l|. In addi- 
tion, the correlations among the mistag parameters and 
those of the At resolution function are reduced. 

The beam-energy substituted mass toes = 

for J/i/jK°) 



V( E bZ m ) 2 ~ (Pes 1 ) 2 ( a11 modes except 
or the difference AE between the candidate center-of- 

mass energy and -E{j™ m (J/tpK° channel) are used to 
determine the composition of our final sample (Fig. [TJl. 
Here, B§™ m and are the beam energy and B 

momentum in the center-of-mass frame. Events with 
toes > 5.2GeV/c 2 (AE < 80MeV) are used so that 
the properties of the background contributions can be 
measured. The more restricted signal region (Table [nj 
contains 7730 CP candidate events that satisfy the 
tagging and vertexing requirements. 

For all modes except r] c K° and J/i/jK® we use simu- 
lated events to estimate the fractions of events that peak 
in the toes signal region due to cross-feed from other de- 
cay modes (peaking background). For the rj c K^ mode 



the cross-feed fraction is determined from a fit to the 
Mkk-k and Toes distributions in data. For the J/ipK® 
decay mode, the composition, effective 77/, and AE dis- 
tribution of the individual background sources are deter- 
mined either from simulation (for B — > J/ipX) or from 
the irig+g— sidebands in data (for fake J/tp — > £ + £~). 

We determine sin2/3 with a simultaneous maximum 
likelihood fit to the At distributions of the tagged Bqp 
and Bfl av samples. The At distributions of the Bqp sam- 
ple are modeled by Eq.^with |A| = 1. Those of the Ba av 
sample evolve according to the known frequency for fla- 
vor oscillation in B° mesons. The observed amplitudes 
for the CP asymmetry in the Bqp sample and for flavor 
oscillation in the Bfl av sample are assumed to be reduced 
by the same factor 1 — 2w due to flavor mistags. The At 
distributions for the signal are convolved with a common 
resolution function, modeled by the sum of three Gaus- 
sians Q. Backgrounds are incorporated with an empir- 
ical description of their At spectra, containing prompt 
and non-prompt components convolved with a resolution 
function |]| distinct from that of the signal. 

There are 65 free parameters in the fit: sin2/3 (1), the 
average mistag fractions w and the differences Aw be- 
tween B° and B° mistag fractions for each tagging cat- 
egory (12), parameters for the signal At resolution (7), 
parameters for CP background time dependence (8), and 
the difference between B° and B° reconstruction and tag- 
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FIG. 1: Distributions for Bcp and Sfl av candidates satisfying 
the tagging and vertexing requirements: a) ttibs for the final 
states J/ipK%, ip{2S)K° s , XdK° s , and r) c K%, b) AE for the 
final state J/ipK° L , c) m ES for J/i>K*°(K*° -» K° s -k°), and d) 
7TIbs for the Bfl a v sample. In each plot, the shaded region is 
the estimated background contribution. 



TABLE II: Number of events iVtag in the signal region after 
tagging and vertexing requirements, signal purity P includ- 
ing the contribution from peaking background, and results of 
fitting for CP asymmetries in the Bcp sample and various 
subsamples. In addition, results on the -Bfl av and charged B 
control samples test that no artificial CP asymetry is found 
where we expect no CP violation (sin2/3 = 0) . Errors are sta- 
tistical only. The signal region is 5.27 < hies < 5.29GeV/c 2 
(\AE\ < lOMeV for J/tpK° L ). 



Sample 


Mag 


P(%) 


sin2/? 


Full CP sample 


7730 


76 


0.722 ± 0.040 


j/y>^,V(2S)^,xcii^ c i^ 


4370 


90 


0.75 ± 0.04 


J/i>K° 


2788 


56 


0.57 ± 0.09 


J/i>K*°(K*° -> K° s ir°) 


572 


68 


0.96 ± 0.32 


1999-2002 data 


3032 


77 


0.74 ± 0.06 


2003-2004 data 


4698 


77 


0.71 ± 0.05 


J/i/>K%, i/>{2S)K%, XciK u s , r, c K» only { Vf = 

TLl.TsUfrsU . _+_-\ niKt nc 


-1) 



J/ipK° s (K° s -» ttV) 653 88 

4>(2S)K° S (K° s — ► TV+TY-) 485 82 

XciK° s 194 81 

r) c K% 287 64 



0.65 ± 0.12 
0.88 ± 0.14 
0.69 ± 0.23 
0.17 ± 0.25 



Lepton category 
Kaon I category 
Kaon II category 
Kaon-Pion category 
Pion category 
Other category 



490 96 

648 93 

1021 89 

769 90 

835 87 

607 88 



0.75 
0.75 
0.77 
0.77 
0.96 
0.23 



± 0.08 
± 0.08 
± 0.09 
± 0.15 
± 0.22 
± 0.51 



-Bfl av sample 


72878 85 


0.021 ± 0.013 


B + sample 


18294 88 


0.003 ± 0.020 



ging efficiencies (7); for -Bfl av background, time depen- 
dence (3), At resolution (3), and mistag fractions (24). 
For the CP modes (except for J/ipK®), the apparent CP 
asymmetry of the non-peaking background in each tag- 
ging category is allowed to float. This asymmetry is pa- 
rameterized so that it does not depend on the value of 
sin2/3. 

We fix t b o = 1.536 ps, Am d = 0.502 ps" 1 [n|, |A| = 1, 
and Ar = 0. The determination of the mistag fractions 
and At resolution function parameters for the signal is 
dominated by the high-statistics -Bfl av sample. Back- 
ground parameters are determined mainly from events 
with m ES < 5.27GeV/c 2 . 

The fit to the Bcp and -Bflav samples yields 

sin2/3 = 0.722 ± 0.040(stat) ± 0.023(syst). 

Figure |21 shows the At distributions and asymmetries in 
yields between B° tags and B° tags for the 77/ = —1 and 
rjf = +1 samples as a function of At, overlaid with the 
projection of the likelihood fit result. 

In a separate fit with only the high purity r/f = — 1 
sample, we obtain |A| = 0.950 ±0.031(stat)±0.013(syst). 
The correlation between the coefficients multiplying the 
sm(ArridAt) and cos(Am,iAi) terms in Eq. ^is —2%. 

The sources of systematic error are summarized in Ta- 




FIG. 2: a) Number of r] f = -1 candidates (J/ifiKl, 
il){2S)K%, XoiK° s , and ri a K%) in the signal region with a B 
tag N B o and with a B tag Ngo , and b) the raw asymmetry 
(N B o — Ng Q )/(N B o + N B o), as functions of At. Figs, c) and 
d) are the corresponding plots for the rjf = +1 mode J/ipK®- 
All plots exclude Other- tagged events. The solid (dashed) 
curves represent the fit projections in At for B° (B°) tags. 
The shaded regions represent the estimated background con- 
tributions. 



bio IIIII These include the uncertainties in the level and 
CP asymmetry of the peaking background, the assumed 
parameterization of the At resolution function, possible 
differences between the -Bfl av and Bcp mistag fractions, 
knowledge of the event-by-event beam spot position, and 
the possible interference between the suppressed b — > ucd 
amplitude with the favored b — > cud amplitude for some 
tag-side B decays |12I . In addition, we include the varia- 
tion due to the assumed values of |A| and AF. We assign 
the change in the measured sin2/3 when we float |A| and 
when we set AT /T = ±0.02, the latter being consider- 
ably larger than recent Standard Model estimates fl3| . 
The total systematic error on sin2/3 (|A|) is 0.023 (0.013). 

The large Bqp sample allows a number of consistency 
checks, including separation of the data by decay mode 
and tagging category, as shown in Table ITU Considering 
statistical errors only, the probability of finding a worse 
agreement in measured sin2/3 values across decay modes 
is 7% and between tagging categories is 86%. The results 
of fits to the control samples of non-CP decay modes 
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TABLE III: Sources of systematic error on sin2/3 and |A|. 


Source 


a(sin2/3) 




CP backgrounds 


0.012 


0.002 


At resolution function 


0.011 


0.003 


J/i/jKl backgrounds 


0.011 


AT / A 

N/A 


Mistag fraction differences 


0.007 


0.001 


Beam spot 


0.007 


0.001 


Am d , tb, AT/r, |A| 


0.005 


0.001 


Tag-side interference 


0.003 


0.012 


MC statistics 


0.003 


0.003 


Total systematic error 


0.023 


0.013 



indicate no statistically significant asymmetry. 

This measurement of sin2/3 supersedes our previous re- 
sult and is consistent with the range implied by other 
measurements and theoretical estimates of the magni- 
tudes of CKM matrix elements in the context of the Stan- 
dard Model 0| . The theoretical uncertainty on the in- 
terpretation of the measurement of sin2/3 in these modes 
is approximately 0.01 [8j- As the current measurement 
is statistics limited, future measurements will add fur- 
ther model-independent constraints on the position of the 
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